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Abstract 

Consistency relations, which relate the squeezed limit of an (N+l)-point correlation function to 
an N-point function, are non-perturbative symmetry statements that hold even if the associated high 
momentum modes are deep in the nonlinear regime and astrophysically complex. Recently, Kehagias 
Riotto and Peloso & Pietroni discovered a consistency relation applicable to large scale structure. We 
show that this can be recast into a simple physical statement in Lagrangian space: that the squeezed 
correlation function (suitably normalized) vanishes. This holds regardless of whether the correlation 
observables are at the same time or not, and regardless of whether multiple-streaming is present. The 
simplicity of this statement suggests that an analytic understanding of large scale structure in the 
nonlinear regime may be particularly promising in Lagrangian space. 
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A Recovering the Eulerian space consistency relation from Lagrangian space — to arbitrary orders 


1 Introduction 


Consistency relations are statements which relate the squeezed limit of an (N+l)-point correlation 
function to an N-point function of cosmological perturbations; i.e., they take the following schematic 
form in momentum space: 


lim 

k^-O 


(7rk0ki0i,,...0k^)^' 

P.(k) 


~ (OkiOk.-.Ok:, 



( 1 . 1 ) 


where vTk represents a squeezed wavemode (long wavelength) of what turns out to be a Goldstone 
boson or pion, Pjr(k) is the power spectrum of the pion (k represents the magnitude of the vector k), 
and O represents observables at high momenta ki,...,kN. The symbol {■■■Y' denotes the connected 
correlation function with the overall delta function removed. Consistency relations can be understood 
as analogues of ‘soft-pion’ theorems in particle physics, which arise generally when a symmetry is 
spontaneously broken/nonlinearly realized. In the case of cosmology, the symmetries in question are 
diffeomorphisms (i.e. coordinate transformations), and consistency relations arise from a particular set 
of residual symmetries of a given gauge where the transformation does not fall off at infinity. The first 
example of a consistency relation was pointed out by Maldacena [3] in the context of a computation 
of the three-point correlation function from inflation. The utility of this as a test of single field/clock 
inflation was emphasized by Creminelli & Zaldarriaga [4]. Recent work pointed out new symmetries 
and therefore further consistency relations [5, 6], indeed an infinite tower of them [7], and explicated 
their non-perturbative nature [8, 9, 10, 11, 12, 13]. 
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These consistency relations are extremely robust: they remain valid when the high momentum 
modes {O in Eq. 1.1) are deep in the nonlinear regime, and even when the observables are astrophys- 
ically complex (such as galaxy density). This point might appear academic when applied to (small) 
perturbations in the early universe, such as are revealed in the cosmic microwave background. When 
applied to large scale structure (LSS) in the late universe, however, the robustness of the consistency 
relations becomes very interesting. It thus came as welcome news when Kehagias & Riotto [1] and 
Peloso & Pietroni [2] (KRPP) pointed out that non-trivial consistency relations exist even if all wave- 
modes (including the squeezed one) are within the Hubble radius, within the Newtonian regime which 
is the realm of LSS (see also [14, 15, 16, 17, 18, 19, 20, 21]). 

The KRPP consistency relation can be stated in the following form: 

^ D(\ • Ic ^ 

( 1 . 2 ) 


lim ' 


(u^(r/)C>ki im) •••Okjv(w))' 


k->-0 


PviK v) 


E (»i) - (w) >“'. 


^0'(,) k2 


where is the j-th component of the peculiar velocity in momentum space, is the velocity 
power spectrum defined by ^k'*(^)) ~ (27r)^(5£)(k — k')(k*k'^7k^)Pi,(k, r/),^ the observables can 

be thought of as mass or galaxy overdensity at different momenta and times, and D and D' represent 
the linear growth factor and its conformal time derivative. The fluctuation variables will in general 
depend on time, although we will often suppress the time dependence to simplify the notation: 

(and its power spectrum) is at conformal time r], Oki is at time r]i, and so on. The times need not be 
equal. The symbol k^ denotes k • k. 

We wish to show that the KRPP consistency relation takes a particularly simple form in Lagrangian 
space: 

(1.3) 

Unless otherwise stated, we use p to denote momentum in Lagrangian space and k to denote momen¬ 
tum in Eulerian space. In other words: 

Ok = yd'xO(x)7‘^-" , Op = y'd3qO(x(q))7P-‘l, (1.4) 

where x and q are the Eulerian space and Lagrangian space coordinates respectively.^ In both cases, 
we rely on context to distinguish between O in Fourier space and O in configuration space. 

Since the velocity v (whose j-th component is v^) is nothing other than the time derivative of the 
displacement A in Lagrangian space, we can also rewrite the Lagrangian space consistency relation 

^This form of the velocity power spectrum assumes no vorticity. This is acceptable since Pu(k) is used only for small k, 

or large scales, where the growing mode initial condition ensures gradient flow. 

^The definitions given apply even in the presence of multiple streaming. See discussion in §2.1. 
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as: 


( 1 - 

where the power spectrum of displacement is defined by (ApA:^,*) = (27r)^(5£)(p — pO(p*pVp^)-Pa(p)- 
It is important to emphasize that the Eulerian space consistency relation (Eq. 1.2) already yields a 
vanishing right hand side if r/i = r] 2 ... = The Lagrangian space consistency relation (Eq. 1.3 
or 1.5), on the other hand, has a vanishing right hand side regardless of what the times r/i,...,r/jv 
happen to be. The consistency relation can also be viewed as a statement about how the squeezed 
correlation function (normalized by the soft power spectrum) scales with the soft momentum: the 
Eulerian space consistency relation states that such a squeezed correlation function goes like (k is 
the soft momentum); the Lagrangian space consistency relation states that there is no p*^ term, and 
at best there is a p^ contribution with e > 0. 

The simplest way to derive Eq. (1.5) is to work out the implications of the KRPP symmetry entirely 
within Lagrangian space. This is done in §2. We perform a perturbative check of this Lagrangian 
space consistency relation using Lagrangian perturbation theory in §3.1. Because the Eulerian space 
and the Lagrangian space relations look so different, as a further check, we show how one can be 
obtained from the other in §3.2. Since observations are performed in Eulerian, not Lagrangian, space, 
the fact that the consistency relation takes a particularly simple form in Lagrangian space is mainly 
of theoretical interest. The simplicity of the Lagrangian space consistency relation should not be 
interpreted as the lack of physical content, however - in the Lagrangian as well as in the Eulerian 
picture, the consistency relation can be viewed as a test of the single-field initial condition and of the 
equivalence principle. Rather, the simplicity suggests that an analytical understanding of nonlinear 
clustering might be most promising in Lagrangian space. This is discussed in §4. 

2 The Lagrangian space consistency relation: derivation 

After a brief review of notation, we derive our main result ~ the Lagrangian space consistency 
relation - using the background wave argument phrased entirely in Lagrangian space. 

2.1 Notation 

We use q to denote the Lagrangian space coordinate of a particle, which coincides with its initial 
position, and x to denote the Eulerian space coordinate which is its position at a later time. To be 
definite, in cases where multiple components are present, the Lagrangian space coordinate q refers 


(Ap(7?)C>pi(??i)...0p^(r/jv))^' ^ 
P^o PA(p,r?) 
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to that of the dark matter particle, which has only gravitational interactions.^ Both coordinates 
are defined in comoving space where the expansion of the nniverse is scaled out. The (dark matter) 
displacement A is the difference: 

x(q,r/) = q +A(q,T/) . (2.1) 


The (dark matter) velocity is given by the conformal time derivative of A at a hxed Lagrangian 
coordinate: 


dA 


( 2 . 2 ) 


The (dark matter) overdensity 5 can be obtained by mass conservation, assuming the initial overdensity 
is negligible: 


1 + (5(x,7/) = |J(q,r/)| ^ 


(2.3) 


with J (q, r]) being the Jacobian relating the volume elements in Eulerian and Lagrangian space: 


J(q, 7]) = det 


dx^{q, Tj) 
'dqi 


(2.4) 


The Jacobian J as a function of q is well-defined even in the presence of multiple-streaming - where 
a single x corresponds to multiple q’s - but Eq. (2.3) requires modification in that case: 


1 + <5(x, 77)= l^(q,f/)l \ (2.5) 

x=q+A(q,r)) 

where the sum is over all q’s that reach the same x. 

Suppose we have some LSS observable O. This could represent many different quantities, such as 
mass overdensity or galaxy number overdensity.^ What we typically observe is O as a function of x 
(and possibly time, which we suppress). Given this function 0(x), one can define unambiguously a 
corresponding function of q: 0(x(q)). In other words, suppose we are interested in the value of O at 
a Lagrangian location q: we can define it by working out the x that q maps to, and then evaluating 
0(x). This procedure is well defined even if multiple q’s map to the same x, which is expected to 
happen for dark matter in the nonlinear regime. 

Some quantities defined in Lagrangian space, on the other hand, might not have an unambiguous 
meaning in Eulerian space. Eor instance, the velocity v given in Eq. (2.2) is defined for a dark matter 

^Our derivation of the Lagrangian space consistency relation would go through even if we chose the Lagrangian coordinate 
to track other constituents of the universe. 

^ Unless otherwise stated, whenever we discuss mass or galaxy density, we mean the mass or galaxy count per unit Eulerian 
space volume. Such a quantity can of course be expressed as a function of either Eulerian space coordinate x or Lagrangian 
space coordinate q. 
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particle labeled by the Lagrangian coordinate q. At an Eulerian position x where multiple Lagrangian 
streams cross, additional inputs are required to define a velocity; a reasonable definition is: 


average ^ v{q), 


x=q+A 


( 2 . 6 ) 


where the sum is over all q’s that map to the same x, and J\f is the number of such q’s. This gives a 
mass weighted velocity. 

It is interesting to contrast the Fourier transform in Lagrangian versus Eulerian space, as described 
by Eq. (1.4). In particular, the Eulerian space Fourier transform can be rewritten as (suppressing time 
dependence): 


0k = y d3x0(x)e*^-^ = J d=^q J(q)0(x(q))e*‘^-('1+^) , 


(2.7) 


where J comes without absolute value; this expression remains valid in the presence of multiple 
streaming. Note how an Eulerian space Fourier transform of 0(x) can be interpreted as a Lagrangian 
space Fourier transform of J(q)0(x(q))e*’^'^. 


2.2 Derivation from the displacement symmetry 

We now deduce our main result, making use of a master formula derived in an earlier paper [21]. 
At the heart of the consistency relation is the existence of a nonlinearly realized symmetry, under 
which some field - the Goldstone boson or pion vr - transforms as vr ^ 7 r +Aiin, 7 r +Ankvr. Here, Aun vr 
is the part of the transformation that is linear in vr, and Ani.vr is the part of the transformation that 
is independent of vr (i.e., nonlinear in tt, though ‘sub-linear’ or ‘inhomogeneous’ would be a better 
description). The fact that Ani.vr 7 ^ 0 is the sign of a nonlinearly realized, or spontaneously broken, 
symmetry. At the same time, there are other fields or observables O that could have their own linear 
and/or nonlinear transformations. The master formula (in momentum space) reads [21]: 

where (...)^ refers to the connected correlation function without removing the overall delta function (as 
opposed to {■■.y which has the delta function removed). Note how it is the nonlinear transformation 
of vr and the linear transformation of O that show up on the left and the right respectively. Note also 
that the 0’s need not even be the same observable. Nor do vr and the 0’s need be at the same time: 
they can be at arbitrary, potentially different, times. The derivation of this master formula made no 
assumption about whether the quantities (or the Fourier transform thereof) are defined in Eulerian 
or Lagrangian space. We are thus free to use it in either. This master relation can be used to derive 
the large scale structure analog of Ward identities or soft-pion theorems in particle physics. 
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As a warm-up, let us first apply this formula to a simple system that involves the dark matter 
only. The dynamics is described by: (1) x = q -|- A as in Eq. (2.1); (2) the dark mater overdensity 6 
determined by the Jacobian as in Eq. (2.5); (3) the displacement A which evolves according to: 


d^A 


drf- 


a' dA 
q a drj 


= , 


(2.9) 


where a is the scale factor, a' is its derivative with respect to conformal time t], is the gravitational 
potential and is the partial derivative with respect to x; lastly (4) the Poisson equation: 




( 2 . 10 ) 


where G is the Newton constant and p is the mean mass density. 
This system has the following symmetry: 


q —> q , A —^ A + n(7/) 






•X, 


( 2 . 11 ) 


where n(7y) is a function of time alone. We will refer to this as the displacement symmetry. Note 
how A shifts by a nonlinear (or sub-linear) amount and can be thought of as our Goldstone boson. 
The same is true for <h. The interesting point is that the mass overdensity 5 does not transform at 
all under this symmetry. Nor are q or t/ transformed. Applying the master formula, choosing the 
observable O = J, we thus hnd: 


1 - 

p^o Pa(p) 


= 0 . 


( 2 . 12 ) 


Here, we have used the fact that the nonlinear transformation of A in Eourier space is A^i.Ap = 
n(ry)(27r)^(5£)(p), where (5_d(p) is the Dirac delta function. We have also removed the overall momentum- 
conserving delta function. The power spectrum of displacement Pa is as dehned in §1. 

Two comments are in order before we proceed to generalize this derivation to more realistic, 
astrophysically complex observables. Eirst, while the Lagrangian coordinate q does not transform 
under the symmetry of interest, the Eulerian coordinate x = q -|- A does, because the displacement 
A shifts. This implies that an observable like 6, when expressed as a function of x, transforms as: 
6 — >■ J-|-Aiin.J with Alin.(5 = A-VJ. Plugging this into the master formula Eq. (2.8), we see that there 
is a non-vanishing right hand side, unlike the situation in Lagrangian space where 6 expressed as a 
function of q does not shift at all. The is the fundamental reason why the KRPP consistency relation 
takes a more complicated form in Eulerian space (Eq. 1.2) than in Lagrangian space (Eq. 2.12). 

Second, the reader might wonder about the validity of our application of the master formula: on 
the one hand, the master relation is phrased in terms of a scalar pion; on the other, our application 
effectively uses the vector displacement A as the pion. The short answer is that the master formula 


7 







is applicable to any field vr that shifts nonlinearly under the symmetry of interest; one can use it for 
each component of A for instance. The long answer is: since A is used in the consistency relation 
only as a soft (long wavelength) mode, one is justified in treating it as a gradient mode (assuming 
growing mode initial condition) with A = Vgvr with vr playing the role of the displacement potential. 

The master formula can then be applied with the displacement potential as the pion. The resulting 
consistency relation can be shown to be equivalent to the one we have derived.^ 

Let us turn to the derivation of a stronger form of the Lagrangian space consistency relation. So 
far, we have focused on a simple system of dark matter particles that interact only gravitationally, as 
embodied in Eqs. (2.9) and (2.10). Let us consider the addition of galaxies into the mix. They have 
their own number overdensity 6g, displacement Ag and velocity Vg = A'g. Their number density is 
not necessarily conserved by evolution, since galaxies can form and merge: 

*^3 “b (1 “b X ' '^g — j (2-13) 

where ' refers to conformal time derivative at a fixed Lagrangian coordinate and Rg is a source term 
that incorporates the formation and merger rates. The equation of motion for the galaxies is: 

A'' + -A'g = -V^^ + Fg, (2.14) 

where encodes additional forces that might act on galaxies, such as gas pressure, dynamical friction 
et cetera. The gravitational potential <1> is determined of course by the Poisson equation (2.10) as 
before. 

The displacement symmetry of Eq. (2.11) can be extended to include also: 

Ag ^ Ag + n(?7), (2.15) 

which also implies Vg —?■ Vg + nb The galaxy overdensity 5g, like its dark matter counterpart, does not 
transform under this symmetry. Eqs. (2.11) and (2.15) represent the displacement symmetry of the 
combined dark-matter-galaxies system, as long as Rg and Fg depend only on (dark matter/galaxy) 
densities and gradients of (dark matter/galaxy) velocities - recall that neither shifts under our sym¬ 
metry. What happens if Rg and/or Fg depends on velocities as opposed to gradients of velocities? 

In that case, shifting velocities by a spatially constant amount would affect the galaxy formation and 
dynamics - this is a violation of the equivalence principle which states that local physical processes 
(such as galaxy formation, mergers and motion) should not be dependent on the absolute state of 

^ There are actually two different nonlinear realized symmetries associated with the displacement potential. One is shifting 
it by a constant or a function of time (but not space). The other is shifting it by a linear gradient, i.e., tt — >■ tt + n ■ q where 
n is the same as that in Eq. (2.11). There are as a result two consistency relations which can be succinctly combined into 
one, Eq. (2.12). See our earlier paper [21] for further discussions. 



motion. Note that a dependence on the dark-matter-galaxy velocity difference v — v^, on the other 
hand, is consistent with the equivalence principle, and the velocity difference is indeed unchanged 
under our symmetry. Thus, as long as the equivalence principle is respected, whether Rg and Fg 
depend on densities, gradients of velocities or velocity differences, the displacement symmetry holds. 
Furthermore, the same statement is expected to be valid in a system with many different species, such 
as baryons, galaxies or even dark matter of different kinds. The argument that leads to Eq. (2.12) 
can be rerun to give the more general Lagrangian space consistency relation: 

( 2 . 

where O is any observable that has no linear shift under the displacement symmetry - this includes 
for instance the densities, displacements and velocities of the galaxies and of dark matter.® Note that 
the d’s need not be the same observables. We have restored the explicit time-dependence of each 
fluctuation variable to emphasize the fact that the times need not be equal. Note also that we have 
chosen the dark matter displacement to be the pion. We could have chosen the galaxy displacement 
instead. Assuming that gravity is the dominant interaction on large scales and adiabatic initial 
conditions, the two displacements are expected to coincide in any case in the soft limit. Furthermore, 
we could have chosen the velocity instead of the displacement as the soft-pion, in which case Eq. (1.3) 
follows. 



3 The Lagrangian space consistency relation: checks 

The above derivation of the Lagrangian space consistency relation is a bit terse, and the form 
the relation takes is surprisingly simple. It is thus worth performing some non-trivial checks of the 
relation. We will first do this using second order Lagrangian perturbation theory (§3.1). Then, in 
§3.2, we demonstrate how the Eulerian space consistency relation can be derived from its counterpart 
in Lagrangian space. 

3.1 Perturbative check 

Let us perform an explicit check of Eq. (2.16) using second-order Lagrangian space perturbation 
theory. For simplicity, we will focus on the case where the only species present is dark matter and the 
observable 0 = 6. We will confine the discussion to the squeezed three-point function; extension to a 

® The reader might wonder: given that the Lagrangian coordinate q does not get transformed at all under the displacement 
symmetry, is there any observable that has a linear shift? The answer is yes. For instance, the combination O = w5 transforms 
to (v -f n')5 giving a shift that is linear in the fluctuation variable 5. 
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general (N+l)-point function is straightforward. Expanding Eq. (2.3) to second order, we have 


1^- .^2 1. 


<5(x(q, rj),rj) = -V^-A + -(Vq • A)^ + - Vq. A*. 


(3.1) 


Expanding out (ip = +(ip^^ H-, Ap = Ap^ + A® H-, and plugging into Eqs. (2.9) and (2.10), 

we have [ 22 ]: 






lD2{r])ip> f d?pAd?VB 


/ 




2D{r]Y p2 
1 f cfpAd^VB 


(2vr) 


-5d(pa + Ps-p) 1- 


/ 


(27r) 


'^d(pa + Pb-p) 1- 


(Pa • Pb? 

pipl 

D 2 {v) , (pa-Pb)^ 
D{r,y 


5pi(^)4s(^)’ 


+ 


p>l 


1 + 


D2ir]) 

D{rif 


(3.2) 


where D is the linear growth factor determining the time-dependence of the first order displacement 
(and density), and D 2 is the second order growth factor determining that of the second order displace¬ 
ment. They satisfy the equations: 

+-L»'- 47rGaV-D = 0 , 
a 

D'i + -D '2 - 4:TTGa^pD2 = -AttGo^pD'^ . (3.3) 

a 

Eor instance, in a flat universe with = I 5 D 2 = —311^/7. Using these expressions, we can work 
out the lowest order contributions to the relevant squeezed bispectrum: 


(^p(^)^Pl (^i)^P2(^2)) 


O(p^) -h 0 (p^Pa(p)) , 

(3.4) 


where we spell out the dependence on the soft momentum p: the 0{yP) piece comes from the first 
term on the right in the first line, and the 0 (P'^T’a(p)) piece comes from the other two terms. We have 
used the fact that 1 — [(pa • Pb)^/p^P^] = O(p^) for Pa + Pb = P- The 0 (p-^(Pa(p)) piece is obviously 
compatible with the Lagrangian space consistency relation; i.e., it gives (Ap( 5 p^( 5 p 2 )^VT’A(p) = 0 in 
the p —)■ 0 limit. The 0{p^) piece does the same, provided the power spectrum Pa(p) is not too 
blue. Parameterizing the power spectrum Ua(p) oc p”~^ in the low momentum limit, the consistency 
relation holds as long as n < 3. Exactly the same condition is needed for the Eulerian space consistency 
relation (see e.g. [ 21 ]).’^ 

^The Eulerian space consistency relation is often given in a form where S is used as the soft mode. One might be tempted 
to do the same for the Lagrangian space consistency relation. However, one can check using perturbation theory that such a 
consistency relation would have required n < 1, a condition considerably stronger than expected. This is related to the fact 
that (5p^^ does not vanish in the p —0 limit, unlike AU^), 
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3.2 Recovering the Eulerian space consistency relation from La- 
grangian space 

The consistency relation takes such a different form in Lagrangian versus Eulerian space that it is 
worth considering how one can be derived from the other. Let us compute the following: 


E — El + Er 


2 

er = -iw Y, 

a=l 


D{r]a) k • ka 
D'{r]) k2 




(3.5) 


The Eulerian space consistency condition is the statement that E = 0. We will content ourselves with 
deriving this - a special case of the more general Eulerian space consistency relation (1.2) - from the 
Lagrangian space consistency relation. 

To relate E to quantities in Lagrangian space, we will slightly abuse our notation. So far, we 
have been using k for the Eulerian space momentum and p for the Lagrangian space momentum. Eor 
instance, i^k^ is defined as 

<5ki = [ d^x5(x) . (3.6) 


Let us rewrite this as 

= I d3qJ(q)5(x(q))e*P-('l+^('l)) 

= J d3qJ(q)5(x(q))e*Pi ‘i + ip’|” J d^q J(q) 5(x(q)) A”^(q) (3.7) 

-\pTPi j d^qJ(q)d(x(q)) A”*(q)A"(q) + ..., 

where J(q) is defined by Eq. (2.4) (with no absolute value). Defining 


d(q) = J(q)(5(x(q)), (3.8) 

we see that the first term on the right of Eq. (3.7) is <5p^, i.e. the Eourier transform of d in Lagrangian 
space. The other terms on the right can likewise be thought of as the Eourier transform of some 
quantity in Lagrangian space. This is why we introduce pi as the momentum label for these Eourier 
components. On the other hand, upon summation, they give the quantity on the left dki=pi which is 
the Eourier transform of density in Eulerian space - this is why we use ki as its momentum label; it 
just happens to take on the numerical value pi which conveniently gives us the appropriate momentum 
label for quantities on the right. It is worth emphasizing that our definitions are general, in that they 
are valid even in the presence of multiple-streaming (see §2.1). 

The expansion in terms of A in Eq. (3.7) is purely formal. In the nonlinear regime, there is no 
sense in which A is small. The expansion provides a convenient way to relate the Eourier transform in 
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Eulerian space to the Fourier transform in Lagrangian space. We will argue E = 0 holds to arbitrary 
order in a power series expansion. 

For the soft mode, we have 


<=P = + •• 


(3.9) 


This is where our abuse of notation is the most egregious: on the left is the velocity Fourier transformed 
in Eulerian space; on the right is the velocity Fourier transformed in Lagrangian space. They agree 
only to lowest order in perturbations. For the soft-mode, ignoring the higher order corrections is 
permissible: the higher order corrections will give higher powers of the soft momentum p compared to 
what is kept in the consistency relation, provided that the soft power spectrum P^(p) or Pa(p) is not 
too blue (see §3.1). Similarly, it can be shown that in the soft limit, there is no need to distinguish 
between in Lagrangian versus Eulerian space.® 

Let us substitute Eq. (3.7) for the hard modes, and Eq. (3.9) for the soft mode, into the expression 
for E in Eq. (3.5). Consider first what contributes to Ei^: 


(^k=p'^kl=pi<5k2=P2) ~('^p^Pl^P2) 


+ 




Pi P2 


('^p^pa(^i)^pi-pa^P2) + 1 2 


d^PA 
(27r)® 

d^PA d^pB 
(27r)3 (27r)3 


/ 


(^^P^Pi- 


Pa^pa 


iVi)dp2- 


A:(4r?2)) 


+ 




d^PA d^pB , ,• j 
(27r)3 (27r)3 \^p^pi-pa- 


PS^PS 


p,A-Ji?i)A^^(r?i)5p,) + (le^2) 


+ 0(A3) + ... 


(3.10) 


where we have largely suppressed the time-dependence to minimize clutter (r/ for the soft mode, and 
rji and r ]2 respectively for the hard modes), except for variables with internal momenta. We emphasize 
that the expansion in A is purely formal, and comes entirely from expanding or The 

first term on the right can be set to zero by virtue of the Lagrangian space consistency condition 
(keeping in mind that this term is divided by T’u(p) as part of the quantity Ei). We will be assuming 
the Lagrangian space consistency relation in its general form (Eq. 1.3): 


lim 

p —>-0 


(Wp(r/) Opi (r?i) ...C»p^(? 7 Ar))^' 


= 0 , 


(3.11) 


^t.(p,^) 

where the observables at hard momenta need not be the same observable. Finally, note that while we 
are interested in the connected part of the correlator on the left hand side of Eq. (3.10), the correlators 


® It is also worth emphasizing that the notion of a well-defined velocity in Eulerian space is valid only when multiple¬ 
streaming is ignored. This is acceptable for the soft-mode. We do not assume single-streaming for the hard modes. 
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on the right hand side are the full correlators, minus the contributions where some proper subset of 
the original hard and soft momenta sum to zero. In particular, the correlators on the right hand side 
of Eq. (3.10) contain both connected and disconnected pieces. 

The second term on the right of Eq. (3.10), formally 0(A), equals 


iV'i 


d^PA 

(27r)3 




+ (^pV-pA)(^PA(m)^P2) + (^P^P2 )(^Pa(^i)^P1-Pa) 


(3.12) 




+ ( 1 ^ 2 ), 


where the connected trispectrum term (vp...)^ (anticipating division by Pvip)) can be set to zero using 
the Lagrangian space consistency relation, the terms involving (up5p^_p^), (vpSp^) and the like have 
one more power of the soft momentum p (and are thus subdominant) compared to terms involving 
(upA™^) which give: 

{2Trf5D{pi + P2 + p)Pv{p,v)ip^^^j^^^^{dpi+p{m)dp2{m)Y' +(1^2). (3.13) 

The third term on the right of Eq. (3.10), formally O(A^), can be treated in a similar way: some 
can be ignored by assuming the Lagrangian space consistency relation, some are subdominant in the 
soft-limit (i.e., they vanish upon division by Pvip) and sending p ^ 0), and the dominant terms are 
those that involve (upA) which give: 


- i27r)^6D (pi + P2 + p)Pv (p, v) 


d^PA 




D'{r]) p 2 


(„))=' + (1 « 2) 


D'irj) p2 


(3.14) 


Thus, combining Eqs. (3.13) and (3.14), El of Eq. (3.5) can be rewritten as: 
i'^l=pi'n) 4i=pi(??i) 42 =p2(%) Y' 


El = lim 

p —>-0 


Pvip,r]) 


= 


+ p"Pi 


DjVi) Pi • P 
D'{ri) p 2 

Pi • P 


{dpiivi)dp2iV2)f - j 


d?PA 

(27r)3 


p^pS* YTI '"IT (4. (')1)V-p.. (’b)A" im)Y 


D'{rj) p 2 


(4i -PA (m ) ^2 iV2) Ap^ im )) 


D'ir]) p2 

+ ( 1 ^ 2 ) + .... 


(3.15) 
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Next, let us rewrite Er using the same strategy: 


Er = - ^ 


D{TJa) P ■ Pa //-„ \ s \ \c 

2 (^ki=pi (^l) Ok2=p2 (^2) ) 


^ D'{r]) p2 


. iD{r]i)p-Pi r, r ( ^^c 
= - TiTTIV —^ (<3pi (^1) <5p2 (r?2)) 


+ 


E'ij]) p2 

d^PA 
(27r)3 


p^p? (4. ('a)4-P^ (’k)a™ ('b))‘ 


(3.16) 


Z)'(r7) p2 


Pi • P, 


+ P^Pr ^ " (<5p,-p,fa)(^p2(r?2)A- fa))- 


+ (le^2) + .... 


D'{rj) p2 

Thus, we see that E = + = 0, at least to the two lowest non-trivial orders in A. The cancelation 

works like this: expanding and as a formal power series in A, a given order for Er is 

canceled by one lower order for Er. It can be shown that this pattern continues to arbitrarily 
high orders. The proof is given in the Appendix. This completes our derivation of the Eulerian 
space consistency relation, embodied in the statement E = 0 (Eq. 3.5), from the Lagrangian space 
consistency relation (Eq. 3.11). 


4 Discussion 


We have shown that the consistency relation takes a particularly simple form in Lagrangian space: 
the squeezed correlation function, suitably normalized, vanishes (Eqs. 1.5): 

(Ap(r/) Op, (rii) ...Op^{riN)Y 


lim 

p —>0 


= 0 . 


(4.1) 


^’A(p,r?) 

where A is the displacement, and O can be many different observables such as mass or galaxy density; 
the quantities can be at different times, and p, Pi, P 2 , ••• label the momenta with p being the soft one.® 
The derivation given in §2.2 is fully non-perturbative and is valid even in the presence of multiple¬ 
streaming. It makes use of a master formula that was derived in an earlier paper [21], which relates an 
(N-l-l)-point function to the linear transformation of an N-point function, for a general nonlinearly- 
realized symmetry (Eq. 2.8). The key realization is that the nonlinearly-realized symmetry of interest 
“ the displacement symmetry - does not require transforming the Lagrangian coordinate q (Eq. 2.11): 


q —>■ q , A —A + 11 ( 77 ) ) —>■ — ( n'' -|-n' 


■ X . 


(4.2) 


where A is the displacement, n is some function of time, is the gravitational potential, a is the 
scale factor and x is the Eulerian coordinate.^® Thus, many observables O such as the mass density 


®See also Eq. 1.3 with velocity v as the soft mode. 

lOff there are multiple species present such as dark matter and galaxies, the same transformation applies to the displacement 
of all species. See Eq. (2.15). 
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or the galaxy density,when expressed as functions of the Lagrangian coordinate q, do not receive 
linear transformations, and so the right hand side of the master formula vanishes.This contrasts 
with what happens when these observables are thought of as functions of the Eulerian coordinate x; 
they receive linear transformations because under the same symmetry 

X —)■ X + n{r}). (4.3) 

It is worth mentioning that the observables O can even be quantities in redshift space. For instance, 
the redshift space mass density (5s(s), where s is the redshift space coordinate, can always be written 
as a function of q, just like what we have done for the real (Eulerian) space mass density. This works 
even in the presence of (real or redshift-space) multiple-streaming. 

It is worth reviewing the assumptions behind the master formula: it assumes (1) single-field/clock 
initial condition - that which follows from single-field inflation - and (2) the adiabatic mode condition 
(first emphasized in [23]) - that the (displacement) mode generated by the symmetry of interest is 
the long wavelength limit of an actual physical mode. In particular, the latter condition requires that 
n(ry) have the same time dependence as the linear growth factor This is the reason why the 

growth factor shows up on the right hand side of the Eulerian space consistency relation (Eq. 1.2). It 
is interesting that because the Lagrangian space consistency relation has a vanishing right hand side, 
the time-dependence of n(r]) has no direct bearing on the form it takes. 

An important point: the fact that the consistency relation can be written in such a simple, even 
trivial, form in Lagrangian space should not be taken to imply the lack of physical meaning. Indeed, 
the consistency relation can be violated if the initial conditions were not of the single-field/clock type. 
Rather, the simplicity suggests that an analytic understanding of nonlinear clustering is perhaps more 
promising in Lagrangian space. This view has a long history, starting from Zeldovich [24] (see also 
[25, 26]). What is interesting is that the consistency relation, by virtue of its being a symmetry 
statement, is non-perturbative, and thus goes beyond perturbative treatments such as the Zeldovich 
approximation. 

Ultimately, observations are performed in Eulerian space, not Lagrangian space. At the nonlinear 
level, the relation between the two descriptions is complex. Our derivation of the consistency relation 
in Eulerian space from its counterpart in Lagrangian space is a case in point (§3.2). It requires a 
formal series expansion in the displacement A. In relating the two descriptions, the expansion in A is 

^^See footnote 4. 

^^Note that even quantities such as A or $ have no linear shift (a shift that is linear in fluctuation variables). More 

complicated observables could have a linear shift; see footnote 6. 

^^In the presence of multiple species, such as dark matter and galaxies, the fact that the displacements for all species have 

to be transformed by the same n(? 7 ) oc D(ri) is a manifestation of the equivalence principle: all species fall at the same rate 
on large scales. 


15 



done in an uneven manner: only phase factors such as are expanded even though other variables, 

such as the density 5, also depend on the displacement. This is not unexpected in relations that are 
purported to be non-perturbative - partial resummation of perturbations is often a useful technique. 
Can our example point to a useful, new resummation scheme? 

A natural question is whether there are relativistic generalizations of statements like Eq. (4.1) - 
consistency relations with a vanishing right hand side. The Lagrangian coordinate (attached to dark 
matter particles) is essentially the freely-falling coordinate. Indeed [27] showed that using the freely- 
falling coordinate, the dilation consistency relation [3] can be rewritten in a similarly simple form 
(see also [28]). Their derivation is perturbative. It should be possible to extend their proof using the 
non-perturbative arguments presented here. More generally, it would be interesting to see if further 
general relativistic consistency relations, such as those found by [7], can also be recast in this fashion. 
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A Recovering the Eulerian space consistency relation 
from Lagrangian space — to arbitrary orders in displace¬ 
ment 


In §3.2, we argue that the Eulerian space consistency relation follows from its Lagrangian space 
counter-part, at least to the two lowest non-trivial orders in a formal expansion in displacement. In 
this Appendix, we show that this works to arbitrary orders. 

We begin by expanding Eq. (3.7) to all orders in A: 


4i=pi = J J(q)(5(x(q))e 

AT . . r 

«=n Pj'i >■■■ ,Pj„ 


ipi-(q+A(q)) 


(A.l) 


A A*1 ... A*" 

°P1-Pm -Pjn'^PU 


where as before, <5(q) = J(q)(5(x(q)). Collecting together the terms of order A"" in the three-point 
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correlator {vp{r])Sp-^^{r|l)5p^{r]2)) in El, we have 


E 

m=0 


ml {n — m)l 


Pi ■■ ■ Pi P2 


■pr 


'Pil:'" ,Pjn 


im) ■ ■ ■ a*p7^ ivi) 


■ '^P2-P1^+1- 


-Pi, 




(A.2) 


where the correlator in the integral is the full correlator containing both connected and disconnected 
pieces, but where no proper subset of the original momenta p,Pi,P 2 sums to zero, since it is the 
connected correlator that appears in El- The correlator can be split into a sum over products of 
connected blocks. Anticipating division by Pv{p), we see that the Lagrangian space consistency 
relation implies that all contributions where the soft velocity U'^(p) is part of a connected correlator 
with two or more other fields will vanish. The remaining terms contain a factor of either 

(^^^(p)4i—•), (^^^(P)5p2-...), (u^(p)Ap.(77i)), {v^{p)Ap-{r]2 )). (A.3) 

The first two types of terms will be suppressed by an additional power of the soft momentum, and 
are subdominant in the squeezed limit. The final set of terms are 


E 


—^ ml (n — m)\ 

m=0 ^ 




•p^ 


'Pjl 5 • ^Pjn 


■ (?7i) • • • Ap7^ {r]i)5p2-pj^^^ 

+ (lo2) 

and they give (relabeling m as m + 1) 

i 


m(u^(r/)A*pi,^ (m))(4i-Pii—-Pi™ im) 

. p , (r,2)Aj77:^ {r^2) ■ ■ ■ K-(m)) 


(A.4) 


.^,.p-piT»(ryi) ^ 

^P _-2 


p 2 D'(r}) m\(n — 1 — m)\ 

^ m=0 ^ ' 


pU...p*™p*^+i 


■p^ 


(^Pl+P-Pil-Pirn (^l) 


'Pjl ’■■■ ’Pin 




Pin 


Pim+1’ 




(A.5) 


+ (1 o2) 

Comparing this to the order terms in the expansion of Eji^ 


n—1 


__9 7 -»//.-\ / ^ 


An—1 


p 2 D'(r}) ^m\(n — l—m)\ 

\ / / V / 


p7...p*-p*™+i 


'P2" 


('^Pl-Pil- 


-Pi 


(^l) 


•A*pb^(r?i)---Aj7^(r,i)4_p^.^^^. 
+ (lf^2) 


' Pil r" ^Pji 

im) ■ ■ ■ K^Av2)) 


(A.6) 


the terms cancel in the p ^ 0 limit to give El + Er = 0 order by order in the formal expansion in A, 

Q.E.D. 


17 








References 


[1] A. Kehagias and A. Riotto, “Symmetries and Consistency Relations in the Large Scale Structure 
of the Universe,” Nucl. Rhys. B 873, 514 (2013) [arXiv: 1302.0130 [astro-ph.CO]]. 

[2] M. Peloso and M. Pietroni, “Galilean invariance and the consistency relation for the nonlinear 
squeezed bispectrum of large scale structure,” JCAP 1305, 031 (2013) [arXiv:1302.0223 [astro- 
ph.CO]]. 

[3] J. M. Maldacena, “Non-Gaussian features of primordial fluctuations in single field inflationary 
models,” JEEP 0305, 013 (2003) [astro-ph/0210603]. 

[4] P. Creminelli and M. Zaldarriaga, “Single field consistency relation for the 3-point function,” 
JCAP 0410, 006 (2004) [astro-ph/0407059]. 

[5] P. Creminelli, J. Noreha and M. Simonovic, “Conformal consistency relations for single-field 
inflation,” JCAP 1207, 052 (2012) [arXiv; 1203.4595 [hep-th]]. 

[6] K. Hinterbichler, L. Hui and J. Khoury, “Conformal Symmetries of Adiabatic Modes in Cosmol¬ 
ogy,” JCAP 1208, 017 (2012) [arXiv:1203.6351 [hep-th]]. 

[7] K. Hinterbichler, L. Hui and J. Khoury, “An Infinite Set of Ward Identities for Adiabatic Modes 
in Cosmology,” JCAP 1401, 039 (2014) [arXiv:1304.5527 [hep-th]]. 

[8] V. Assassi, D. Baumann and D. Green, “On Soft Limits of Inflationary Gorrelation Functions,” 
JCAP 1211, 047 (2012) [arXiv:1204.4207 [hep-th]]. 

[9] V. Assassi, D. Baumann and D. Green, “Symmetries and Loops in Inflation,” JHEP 1302, 151 
(2013) [arXiv:1210.7792 [hep-th]]. 

[10] A. Kehagias and A. Riotto, “Operator Product Expansion of Inflationary Gorrelators and Gon- 
formal Symmetry of de Sitter,” Nucl. Phys. B 864, 492 (2012) [arXiv:1205.1523 [hep-th]]. 

[11] W. D. Goldberger, L. Hui and A. Nicolis, “One-particle-irreducible consistency relations for 
cosmological perturbations,” Phys. Rev. D 87, 103520 (2013) [arXiv: 1303.1193 [hep-th]]. 

[12] G. L. Pimentel, “Inflationary Consistency Conditions from a Wavefunctional Perspective,” JHEP 
1402, 124 (2014) [arXiv;1309.1793 [hep-th]]. 

[13] L. Berezhiani and J. Khoury, “Slavnov-Taylor Identities for Primordial Perturbations,” JCAP 
1402, 003 (2014) [arXiv; 1309.4461 [hep-th]]. 

[14] P. Creminelli, J. Noreha, M. Simonovic and F. Vernizzi, “Single-Field Consistency Relations of 
Large Scale Structure,” JCAP 1312, 025 (2013) [arXiv; 1309.3557 [astro-ph.CO]]. 

[15] M. Peloso and M. Pietroni, “Ward identities and consistency relations for the large scale structure 
with multiple species,” JCAP 1404, 011 (2014) [arXiv; 1310.7915 [astro-ph.CO]]. 


18 



[16] A. Kehagias, J. Norefia, H. Perrier and A. Riotto, “Consequences of Symmetries and Consistency 
Relations in the Large-Scale Structure of the Universe for Non-local bias and Modified Gravity,” 
Nucl. Phys. B 883, 83 (2014) [arXiv;1311.0786 [astro-ph.CO]]. 

[17] P. Berger, A. Kehagias and A. Riotto, “Testing the Origin of Cosmological Magnetic 
Fields through the Large-Scale Structure Consistency Relations,” JCAP 1405, 025 (2014) 
[arXiv;1402.1044 [astro-ph.CO]]. 

[18] P. Creminelli, J. Gleyzes, M. Simonovic and F. Vernizzi, “Single-Field Consistency Relations 
of Large Scale Structure. Part IL Resummation and Redshift Space,” JCAP 1402, 051 (2014) 
[arXiv;1311.0290 [astro-ph.CO]]. 

[19] P. Valageas, “Consistency relations of large-scale structures,” arXiv:1311.1236 [astro-ph.CO]. 

[20] P. Creminelli, J. Gleyzes, L. Hui, M. Simonovic and F. Vernizzi, “Single-Field Consistency Re¬ 
lations of Large Scale Structure. Part III: Test of the Equivalence Principle,” arXiv; 1312.6074 
[astro-ph.CO]. 

[21] B. Horn, L. Hui and X. Xiao, “Soft-Pion Theorems for Large Scale Structure,” arXiv; 1406.0842 
[hep-th]. 

[22] F. Bernardeau, S. Colombi, E. Gaztanaga and R. Scoccimarro, “Large scale structure of the 
universe and cosmological perturbation theory,” Phys. Rept. 367, 1 (2002) [astro-ph/0112551]. 

[23] S. Weinberg, “Adiabatic modes in cosmology,” Phys. Rev. D 67, 123504 (2003) [astro- 
ph/0302326]. 

[24] Y. B. Zeldovich, “Gravitational instability: An Approximate theory for large density perturba¬ 
tions,” Astron. Astrophys. 5, 84 (1970). 

[25] R. A. Porto, L. Senatore and M. Zaldarriaga, “The Lagrangian-space Effective Field Theory of 
Large Scale Structures,” JGAP 1405, 022 (2014) [arXiv:1311.2168 [astro-ph.CO]]. 

[26] L. Senatore and M. Zaldarriaga, “The IR-resummed Effective Eield Theory of Large Scale Struc¬ 
tures,” arXiv: 1404.5954 [astro-ph.CO]. 

[27] E. Pajer, F. Schmidt and M. Zaldarriaga, “The Observed Squeezed Limit of Cosmological Three- 
Point Functions,” Phys. Rev. D 88, no. 8, 083502 (2013) [arXiv: 1305.0824 [astro-ph.CO]]. 

[28] T. Tanaka and Y. Urakawa, “Dominance of gauge artifact in the consistency relation for the 
primordial bispectrum,” JCAP 1105, 014 (2011) [arXiv:1103.1251 [astro-ph.CO]]. 


19 



